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Abstract 

The recently developed theory of quasi-Lie schemes is studied and 
applied to investigate several equations of Emden type and a scheme 
to deal "with them and some of their generalisations is given. As a first 
result we obtain t-dependent constants of the motion for particular 
instances of Emden equations by means of some of their particular 
solutions. Previously known results are recovered from this ne"w per- 
spective. Some t-dependent constants of the motion for equations of 
Emden type satisfying certain conditions are recovered. Finally ne"w 
exact particular solutions are given for certain cases of Emden equa- 
tions. 

1 Introduction. 

Systems of nonautonomous first-order differential equations appear broadly 
in Mathematics, Physics, Chemistry and Engineering. Therefore methods to 
solve these systems and analyse their properties are specially interesting be- 
cause they allo"w us to understand many important problems in these various 
fields. 

As a first insight into this topic the theory of Lie systems can be con- 
sidered [U [21 E]. Many applications have recently been studied through this 
theory [H [5] . Nevertheless there are many differential equations which cannot 
be studied by means of the theory of Lie systems [HI C] . To treat more general 
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differential equations a new theory has been recently developed. This theory 
is based on the so-called quasi-Lie schemes [B]. These schemes generalise the 
concept of Lie system and sometimes they allow us to transform a system of 
differential equations, generally a non-Lie system, into a Lie system. Once 
the final Lie system is studied and its properties are found, the theory of 
quasi-Lie schemes provides constants of the motion, t-dependent superposi- 
tion rules or even solutions for the initial system of differential equations. 

Quasi-Lie schemes have been useful to deal with many systems of differ- 
ential equations. They have been applied to study nonlinear oscillators 
dissipative Ermakov systems [H [8], etc. In this paper we apply quasi-Lie 
schemes to investigate the properties of Emden-Fowler equations. The lit- 
erature about these equations is very large [n]-[IS] and applications of these 
equations can be found, for example, in Mathematical Physics, Theoretical 
Physics, Astronomy, Astrophysics and Chemical Physics. We only consider 
a small sample and, for instance, there are about 140 references in the review 
by Wong in 1977 [HI . 

Our aim in this paper is to investigate the properties of Emden-Fowler 
equations from the point of view of the theory of quasi-Lie schemes and 
quasi-Lie systems. We firstly show that the knowledge of certain partic- 
ular solutions allows us to transform a given Emden-Fowler equation into 
a Lie system by means of a quasi-Lie scheme and determine a constant of 
the motion through this particular solution. We also study the generalised 
Emden-Fowler equation to recover from our point of view the origin of the 
Kummer-Liouville transformation. Next, use is made of the transformation 
properties of our quasi-Lie scheme in order to obtain some t-dependent con- 
stants of the motion for families of Emden type equations satisfying some 
conditions. Finally, we derive a family of exact particular solutions for a 
particular Emden-Fowler equation by means of a certain kind of t-dependent 
superposition rule. 

In this paper we start in Section 2 with a report on the theory of Lie 
systems and Lie schemes in order to review in Section 3 some previous re- 
sults about quasi-Lie schemes and Emden equations and provide some new 
details. Section 4 is devoted to show that certain particular solutions of Em- 
den equations enable us to obtain t-dependent constants of the motion. The 
latter result is applied in Section 5 to study some particular cases of Em- 
den equations and build up some particular instances of such equations for 
which we can obtain a t-dependent constant of the motion by means of our 
method. A generalised Emden-Fowler equation is studied in Section 6 and 
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we get the Kummer-Liouville transformation from our framework. In Section 
7 we use the transformation properties of our scheme to obtain t-dependent 
constants of the motion for certain equations of Emden type satisfying some 
integrabihty conditions. Finally, in Section 8 we apply some results obtained 
along the paper to analyse certain Emden- Fowler equations. As a result, we 
recover some results about these equations and we find a t-dependent partial 
superposition rule. Such a superposition rule is used next to obtain a family 
of solutions for an important particular Emden-Fowler equation. 



2 Review on Lie and quasi-Lie systems. 

In this Section we report some previous results on the theory of Lie 
systems, t-dependent vector fields and quasi-Lie systems and schemes. Fur- 
thermore we define the extended group of transformations in order to improve 
the methods of the theory of quasi-Lie schemes. Most results and mathemat- 
ical objects used here are only locally defined, but for the sake of simplicity 
we drop this kind of technical detail for the time being. For a full description 
of the basic details see [21 E] • 

A nonautonomous system of first-order ordinary differential equations in 
a manifold is represented by a t-dependent vector field X = X(t,x) on 
such a manifold. The system of differential equations associated with the 
t-dependent vector field X{t, x) is written in local coordinates as 

— = X''{t,x), i = 1, . . . ,n = dimN, 

where X{t, x) = Xlj'Li ^*(^) x)d/dx\ The conditions for this system ensuring 
that it admits a superposition rule, i.e. there exists an open U C A^'"+^ and 
a map $ : f/ C A^("^+^) — >• A^ such that its general solution can be written as 

X(t) = . . . ,X(rn)(t);ki, . . . , kn), 

where \ a = 1, . . . ,m} is a family of particular solutions and k = 

{ki, . . . , kn) is a set of n arbitrary constants such that 

ki,...,kn) C U, 

were studied by S. Lie The necessary and sufficient condition is that 
the associated t-dependent vector field, X(t), can be written as a linear 
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combination 

r 

X(t) = J]&„(t)X(„), (1) 

Q=l 

where the vector fields \ a = 1, . . . ,r} are hnearly independent vector 

fields, i.e. if Ai, . . . , are real constants such that J2a=i ^aX(a) = 0, then 
Ai = • • • = Ar = 0, generating the so-called Vessiot-Guldberg Lie algebra V 
of vector fields. The latter also means that there exist real numbers, 
such that 

r 

[X(„), X(^)] = ^ Ca/3^X(-^) , a, P = 1, . . . ,r. 

7=1 

Any Lie system described by a t-dependent vector field on a manifold X, 
like ([T]) , where the vector fields are complete and close on a finite-dimensional 
Lie algebra determines a free, but maybe a discrete set of points, left action 
$:GxX— >Xofa Lie group G with Lie algebra g ~ \^ on the manifold 
X describing the vector fields in V as fundamental vector fields of such an 
action. 

A solution of this system is represented by a curve s 7(5) in X (integral 
curve) the tangent vector of which 7 at t satisfies 

7(t)=X(t,7(t)). (2) 

It is well-known that, at least for the smooth X with which we work, for 
each xo there is a unique maximal solution 7^ (^) of system ([2]) with the 
initial value xq, i.e. satisfying 7^(0) = xq. The collection of all maximal 
solutions of the system ([2]) gives rise to a (local) generalised fiow g-^ on X. 
By a generalised flow g on N we understand a smooth t-dependent family 
gt of local diffeomorphisms on X, gt{x) = g(t,x), such that go = 16.^. The 
generalised fiow g-^ induced by the t-dependent vector field X is defined by 

^^^(t,a;o)=7x''W- (3) 
Note that for g = g^ equation ([3]) can be formally rewritten in the form 

X{t)=X{t,x)=gtog-K (4) 

We observe that equation (jlj) in fact defines a one-to-one correspondence 
between generalised fiows and t-dependent vector fields. Any two generalised 
fiows g and h can be composed: by definition {goh)t = gtohf. As generalised 
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flows correspond to t-dependent vector fields, this gives rise to an action of a 
generalised fiow hon a t-dependent vector field X, giving rise to h^X defined 
by the equation 

= hog^. (5) 
A more explicit form of this action is 

{h^X)t = htoh-' + {htUX{t)), (6) 

where (/it)* is the standard action of diffeomorphisms on vector fields. These 
results can be summarised by means of the following theorem. 

Theorem 1. The equation ^ defines a natural action of generalised flows 
on t-dependent vector fields. This action is a group action in the sense that 

{goh)^X = g^{h^X). 

The integral curves ofh^X are of the form ht{'~f{t)) for'~f{t) being an arbitrary 
integral curve for X . 

We must notice that ([6]) and the theorem above still works even if h is 
a one-parameter set of diffeomorphisms hf : N ^ N with Hq ^ Id and we 
can define the t-dependent vector field h^X as the one with integral curves 
ht{'~f{t)), where 7(t) is any integral curve for X. 

Let us state the fundamental concepts of the theory of quasi-Lie schemes. 

Definition 1. Let W and V be nonnuU finite-dimensional real vector spaces 
of vector fields on a manifold A^. We say that they form a quasi-Lie scheme 
S{W, V) if they fulfil the conditions below: 

1. W is a. vector subspace of V. 

2. is a Lie algebra of vector fields, i.e. [W, W] C W. 

3. W normalises V, i.e. [W, V] C V. 

If is a Lie algebra of vector fields V, we call the quasi-Lie scheme S(y, V) 
simply a Lie scheme S{V). 

There is the largest Lie subalgebra we can use as W - the normalizer of 
V in V. Sometimes, however, it is useful to consider smaller Lie subalgebras 
W. 
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We say that a t-dependent vector field X is in a quasi-Lie scheme S{W, V) 
and write X G S{W, V) if X belongs to V on its domain, i.e. X(t) G V. 

Now, given a quasi-Lie scheme S{W, V) which we call sometimes simply 
a scheme, we may consider the group, G(W), of generalised fiows associated 
with W. 

Definition 2. We call the group of the scheme S{W, V) the group Q{W) of 
generalised fiows corresponding to the t-dependent vector fields with values 



Given a scheme, S{W, V), with W a Lie algebra of complete vector fields 
we can associate W with a left action ^ : G x N —>■ N, with TeG ~ W, 
describing their elements as fundamental vector fields. Then the generalised 
fiows of Q(W) are those of the form gt{x) = ^{g{t),x) = <l>g(f)(x) with g{t) a 
curve in G with g{0) = e; see [3]. 

Proposition 1. Given a scheme S{W,V), a t-dependent vector field X G 
S{W, V) and a generalised fiow g G GiW), we get that g^X G S{W, V). 

We look for a set of t-dependent transformations containing Q{W) satis- 
fying the latter proposition. 

Lemma 1. Consider a scheme S{W,V) with W a Lie algebra of complete 
vector fields. Given an element g G exp(0) and a vector field X G S{W, V), 
then <!?g,X G S{W,V). 

Proof. As (7 G exp(g), there exists an element a G g such that g = exp(a). 
Consider the curve /i : t G [0, 1] — >■ exp(ta) G G. By means of the action 
^ : G X N N associated with the Lie algebra of vector fields W of the 
scheme S(W,V), the curve h(t) induces the generalised fiow : {t,x) G 
Rx N ^ h({x) = $(exp(t ) G N for the vector field 



in W. 



d 
dt 



i=0 




d 
dt 



t=0 



$(exp(t a), x) 



and Y G W. For each t define the vector field Zj: 



h^X to get 




6 



If we call ' — hJ^{[Y, X]) and apply the last formula to [Y, X], we get 

ihl{[Y,X]))^ = [Y,Xl + ^Z^^\x)ds' 

^[Y,Xl+ f\hJ,,[Y,[Y,X]]Us'. 
Jo 

Defining Z^''^ in an analogous way and applying all these results to the initial 
formula for hJ^X we obtain 

(hlX), ^ X, + [Y,Xlt + ^[Y, [Y,X]lt' + [Y, [Y,X]]lt' + ■■■. 

By means of the properties of the scheme we obtain that each term belongs 
to the scheme, i.e. 

[Y,[Y,...,[Y,X]...]]eS{W,V), 

and therefore 

^g,X^h^,X e S{W,V). 

□ 

Proposition 2. Consider a scheme S{W, V) with W a Lie algebra of com- 
plete vector fields. Given a curve g{t) C exp(g) and a t-dependent vector 
field X{t) E S{W, V), then g^X e S{W, V). 

Proof. It has been shown that 

{9-kX)t^9togr^ + 9t*{^), 

but gt o g^^ E W E S{W, V) and by means of the Lemma 1 we get that 
gt*X E S{W, V) for each t. Hence we have g^^X E S{W, V). □ 

Definition 3. Given a scheme S{W, V) we call the extended group of the 
scheme, Ext(W^), the set of t-dependent transformations ^g{t) induced by 
curves g{t) C exp(0) and the action $ associated with the Lie algebra of 
complete vector fields W. 

Prom the last definition we can state the definition of quasi-Lie system 
with respect to a scheme. 

Definition 4. Given a quasi-Lie scheme S{W,V) and a t-dependent vector 

field X E S{W,V), we say that X is a quasi-Lie system with respect to 
S{W, V) if there exists a t-dcpcndent transformation g E Ext{W) and a Lie 
algebra of vector fields Vq gV such that 

g^X E S{Vo). 
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3 The Emden equation. 

In this Section we approach from the perspective of the theory of quasi-Lie 
schemes the so-called Emden equations of the form 

X = a{t)x + n^l. (7) 

These equations can be associated with the system of first-order differential 
equations 

i) = a{t)v + b{t)x''. 

This system was already studied in [6J from the point of view of the theory 
of quasi-Lie schemes. We summarise next the results of that paper and we 
use them to obtain new properties: t-dependent constant of the motion by 
means of particular solutions, reducible particular cases of Emden equations 
etc. 

Consider the real vector space, VEmd, spanned by the vector fields 

Xi = xdy, X2 = x^dy, Xs = vdr,, X4 = vdv, X5 = xdx- 

The t-dependent vector field determining the dynamics of system ([8]) can be 
written linear combination 

X{t) = a{t)Xi + X^ + h{t)X2. 

Moreover the linear space VFEmd C VEmd spanned by the complete vector 
fields, 

Fi = X4 = vdy, ¥2= Xi = xd^, F3 = X5 = xd^, 

is a three-dimensional real Lie algebra of vector fields with respect to the 
ordinary Lie Bracket because these vector fields satisfy the relations 

[Y,,Y2] = -Y2, [YuYs] =0, [^2,^-3] = -^2. 

Also [WEmd, Vsmd] C Vsmd bccause 

[Fi, X2] = -X2, [Y^, X3] = X3, [F2, X2] = 0, 

[F2,X3]=X5-X4, [Y3,X2]=nX2, [1^3,^3] =-X3. 
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So we get a quasi-Lie scheme SiWEmd, ^md) which can be used to treat the 
Emden equations ([HD • This suggests that we perform the t-dependent change 
of variables associated with this quasi-Lie scheme, namely, 



7(t),/3(t)>0,Vt, (9) 



V = f3{t)v' + a{t)x' , 
which transforms the original system into 

( dx^^fa{ti_j{ty\ m , 

dt VtW lit)) ' 

- 'ait) - 44 - + 44 (ait) - 44 - 44 + 44 ) X' (10) 



dt \" j{t) m j m \ ' ' lit) ait) Jit) 
I mrit) ,n 



The key point of our method is choosing appropriate functions, a, /3 and 
7, in such a way that the system of differential equations ffTOl) becomes a Lie 
system. 

A possible way for the system ( |TOl) to be a Lie system is to choose functions 
a, P and 7 such that the latter system is determined by a t-dependent vector 
field Xit) = fit)X, where X is a true vector field and fit) is no vanishing 
function in the intervel of t under study. As is shown in next Section, this 
cannot always be done and some conditions must be imposed on the initial 
i-dependent functions, a,P and 7, assuring such a transformation to exist. 
These restrictions lead to integrability conditions. 

Suppose for the time being that this is the case. Therefore the system 
is 

-TT = fit) icnx' + Cuv') , 

d.' (11) 

= fit)ic2iv' + C,X' + C22X"') 

and it is determined by the t-dependent vector field 

Xit) = fit)X, 

with 

X = (ciix' + Ci2v')dr,' + (c22a;'" + c^x' + C2iv')d^'. 
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Under the t-reparametrisation, 



Tit) = J /(tOdf, 

system flTT]) is autonomous. The new autonomous system of differential equa- 
tions is determined by the vector field X on TM and therefore there exists a 
first integral. This can be obtained by means of the method of characteristics, 
which provides the characteristic curves where the first-integrals for such a 
vector field X are constant. These characteristic curves are determined by 

dx' dv' 



CnX' + CuV' C2lV' + CxX' + C22X'"- 

which can be written as 

{c2iv' + Cxx' + C22a;'")dx' - {cux + Ci2v')dv' = 0. (12) 
This expression can be straightforwardly integrated if 

dv'{c2iv' + c^x' + C22a;'") = -d^'^Cnx' + Ci2v') =^ C21 = -Cu- (13) 
Under this condition we obtain the first integral for f|T2|) . namely 



I = -Cl2— + C^— + C2lV' X' + C22-^^-^- (14) 

Finally, if we write the latter expression in terms of the initial variables x, v 
and t, we get a constant of the motion for the initial differential equation. 

When we do not want to impose condition f[T^ . we could also integrate 
equation (fT2|) by means of an integrating factor, namely, we look for a func- 
tion, /i(x',f'), such that 

dv' {fi{c2iv' + c^x' + C22a;'")) = dxi{-^[cx\x' + Cv2v')). 

Thus the integrating factor satisfies the partial differential equation 

^(c2if' + c^x' + C22a;'") + ^(ciia;' + Cx2v') = -Ai(cii + C21). 

If Cii + C21 = 0, the integral factor can be fixed to be /i = 1 and we get the 
latter first integral f|T^ . On the other hand, if Cn + C21 7^ 0, we can still look 
for a solution for the partial differential equation for fi and obtain a new first 
integral. 
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4 t-dependent constants of the motion and 
particular solutions. 

Our aim is to show that the knowledge of a certain particular solution of 
the Emden equation satisfying a certain condition allows us to transform it 
into a Lie system and to obtain a t-dependent constant of the motion. Even 
if some results related to this t-dependent constants of the motion can be 
found in the literature, here we recover them from a new point of view. 

If we restrict ourselves to the case a{t) =0 in the system of differential 
equation (ITUl) . it reduces to 

^ = _M ' , M ' 

dt r^'^ Pit) r Pit) 



In order to transform the original Emden-Fowler differential equation 
into a Lie system by means of our quasi-Lie scheme, we try writing the 
transformed differential equation in the form 

-TT = fit) icux' + cuv') , 
dt 

where the Cij are constants. This system of differential equations can be 
reduced to an autonomous one because under the t-dependent change of 
variables 

r(t)= r/(t')dt' 



the latter differential equation becomes 

— = Cux' + Ci2f ', 

— = C22X + C21V . 
dr 

In order for system (fT5l) to be like system (fT6l) we look for functions a, P and 
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(18) 



7 satisfying the conditions, 

The conditions of the first hne lead to 

Pit) = -^7W, (19) 
Cu 

and using this equation in the last relation we obtain 

m = ^ - (20) 

C21 C21 7(t) 

On the other hand, from the three first relations in f|T8l) we get 

C22C12 71^) 

The equality of the right hand sides of (!20l) and (pTI) leads to the following 
equation for the function 7: 

J = a(t)7 + £li^6(t)7-. 
C22C12 

Suppose that we make the choice, with C21 = — cn as indicated in (fT3l) . 

C22 = -1, cii = 1, C21 = -1, C12 = 1 (22) 

and thus (ciiC22)/(c2iCi2) = 1. Therefore we find that 7 must be a solution 
of the initial equation (JTj). In other words, if we suppose that a particular 
solution Xp{t) of the Emden equation is known, we can choose •yit) = Xp{t). 
Then, according to the expression f|T9|) and our previous choice (|22|) . the 
corresponding function [3 turns out to be 

(3{t) = -xp{t). 

Finally, in view of conditions f|T8|) . we get that 



Clll{t) C22/3(t) 

12 



and taking into account our choice (12^ and 7(t) = Xp{t), we obtain the 
condition satisfied by the particular solution: 

xT\t)=xlit). (23) 

The system of differential equations ( fT6l) for such a choice ( l22i) of the con- 
stants {cij = 1,2} is the equation for the integrals curves for the t- 
dependent vector field 

X{t) = fit) {{x' + v') d,, - [v' + a;'") d^^) . 

The method of the characteristics can be used to find the following first- 
integral for this vector field and, in view of f|T^ . we get 



l{x\ v') = + + a^V, n i {-1, 1}, 

/(x', v') = log x' + -f + x'v', n = — 1, 

and, if we express this integral of motion in terms of the initial variables and 
t, we obtain a, as far as we know, new t-dependent constant of the motion 
for the initial Emden equation 

X""*"-^ V'^ XV J r . n 

{n + l)x;+\t) 2xl{t) Xp{t)xp{t) 

/ \ 2 V^'*/ 

r/ ^ I X \ V XV 

I{t,x,v) = \og\—-\+—^ , n = -l. 

\Xp{t) J 2x^{t) Xp[t)xp[t) 

So the knowledge of a particular solution for the Emden equation holding fl23!) 
enables us first to obtain a constant of the motion and then to reduce the 
initial Emden equation into a Lie system. Thus all Emden equations with 
such particular solutions are quasi-Lie systems with respect to the above 
mentioned scheme, the applicability of which depends on a prior knowledge 
of such a particular solution. 



5 Applications of particular solutions to Em- 
den equations. 

This Section is devoted to illustrate the theory with some particular in- 
stances of Emden equations for which one is able to find a particular solution 
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satisfying an integrability condition in an easy way and use is made of such 
a solution in order to obtain the corresponding t-dependent constant of the 
motion. In some cases the so obtained constant can be found in the hterature, 
but they are found here from a new systematic procedure. 

We start with a particular case of the Lane-Emden equation 

2 

x = -jx-x^. (25) 

The most general Lane-Emden equation is generally written as 

2 

X = + f{^) 

and the example here considered corresponds to f{x) = —a;", n ^ 1, which 
is one of the most interesting cases, together with that of f{x) = —e~^^. 
Equation (l25l) appears in the study of the thermal behavior of a spherical 
cloud of gas |tl8j and also in astrophysical applications. A particular solution 
for (l25l) satisfying (l23l) is Xp {t) = (2t)-^/2_ jf substitute this expression 
for Xp{t) and the corresponding one for Xp{t) into the t-dependent constant 
of the motion fl2^ . we get that 

r{t,x,v) = + U^v^ + Afxv (26) 

3 

is a t-dependent constant of the motion proportional to and also pro- 
portional to the t-dependent constants of the motion found in [6l HH [20] . 

We study from this new perspective other Emden equations investigated 
in [T6]. Consider the particular instance 

5.2 

X = X — X . 

t + K 

A particular solution for this Emden equation satisfying (l23i) is 

4 

Xp{t) 



(t + Kf 

In this t-dependent constant of the motion is 

I\t,x,v) = lx^{t + Kf + ^v\t + Kf + 2xv{t + Kf 
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which is proportional to the one found by Leach in 
Now another Emden equation found in 



3.9 

2(t + K) 

admits the particular solution 

- 1 

which satisfies ( 123|) . The corresponding t-dependent constant of the motion 
is given by 

/'(t, X, v) = {K + t)3/2(10(if + t)f 2 + 5vx + 2{K + 

which is proportional to that given in [16] . 

As another example of Emden equation we can consider 

5 . 7 

rr* /y» 

3(t + K) 

which admits as a particular solution 

1 



Xp{t) 



3V3(t + i^)l/3' 

which satisfies (1231) and leads to the t-dependent constant of the motion 

I'{t, X, v) = {K + t)^/%12{K + t)v'^ + 8vx + 3x\K + t)). 

Finally we apply our development to obtain a t-dependent constant of 
the motion for the Emden equation 

with 

n + 3 

We can find a particular solution of the form 
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In order Xp{t) to be a particular solution we must hold the following relation 



{Ki + K3ty+^ {Ki + Ksty+^ {Ki + K^ty- 

and thus 

z/ + 2 = nz/ and v{v + 1)KIK2 = vK2K^ - Kl^. 
From these equations we get 

2 1 2^ 

V = K = 

n — 1 ' {n + Sy 

Under these conditions it can be easily verified that Xp(t) = Xp'^^it). Thus a 
i-dependent constant of the motion is 



n + 1 2 J 



+ {Ki + ir3t)("+=^)/("-^)^^, (28) 

which can be found also in ^6]. 

Another advantage of our method is that it allows us to obtain Emden 
equations admitting a previously fixed t-dependent constant of the motion. 

Suppose we want to construct an Emden equation with a previously cho- 
sen particular solution, Xp(t), satisfying x^it) = Xp'^^{t) for certain n e 
M — {1, — 1}. We can integrate this equation to get all possible particular 
solutions which can be used by means of our method, i.e. 

, 2_ 

1 — n ^ 

We consider functions a{t) and b{t) such that 

Xp = a{t)xp + &(t)Xp. 
For the sake of simplicity we can suppose that h{t) = — 1. Then we get 

P ' p 



a(t) = 

16 



Xp 



If we substitute in this expression for a{t) the chosen particular solution, we 
obtain 

(+\ - 3 + n 

- 2{K + '-^ty 

which leads to an Emden equation equivalent to (1271) and the t-dependent 
constant of the motion for this equation is again (l28l) . In this way we recover 
the cases studied in this Section. 



6 The Kummer-Liouville transformation for 
a general Emden-Fowler equation. 

The general form of the Emden-Fowler equation considered nowadays is 

x + p{t)x + q{t)x = r{t)x''. (29) 

This generalisation arises naturally as a consequence of our scheme. This 
latter second-order differential equation is associated with the system of first- 
order differential equations 



X = V, 

V = —p{t)v — q(t)x + r(t)x"', 



(30) 



which is the system for the determination of the integral curves for the t- 
dependent vector field 

X{t) = -p{t)X, - q{t)X, + r(t)X2 + X3. 

This t-dependent vector field is a more general case than the one studied in 
previous Sections. Under the set of transformations ^ the initial system 
(|30|) becomes the new system 



ck/ 



lit) lit) 



lit) 



-Pit) - 44 - §4 + #1 i-vit) - 4$ - 44+ 



lit) m 



lit) a{t) 



lit) 



a{t) 
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If we choose a = 7, the system reduces to 
dt ~ j{t) ' 



'Pit) 



lit) Pit) 



v + 



Pit) 



-pit)-M-,it)M\x' 



at) 



at) 



I rit)rit) ,n 

^ Pit) • 



When the function 7(t) is chosen to be such that 7 = —qit)'y — p(t)7, i.e. 7 
is a solution of the associated linear equation, we obtain 



r dx' 
dt 



lit) ' 



lit) Pit) 



(31) 



Finally, if the function Pit) is such that 

at) Pit) 



we obtain 



-Pit) - 



dx^ 

Id 



lit) Pit) 



Pit), 



0, 



lit) 



dz/ ^ r{t)rit) ,n 
I dt f3{t) ' 

which is related to the second-order differential equation 



(32) 



d^x' , ,7 



with 



Tit) 



P\t) 
'Pit' 



-X 



In 



lit') 



dt'. 



The new form of the differential equation is called the canonical form of the 
generahsed Emden-Fowler equation. 

This fact is obtained in the previous literature by means of an appropriate 
Kummer-Liouville transformation, but here we obtain it as a straightforward 
application of the properties of transformation of quasi-Lie schemes thereby 
underscoring the theoretical explanation of such a Kummer-Liouville trans- 
formation. 
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7 Constants of the motion for systems of Emden- 
Fowler equations. 

In this Section we show that under certain assumptions for the t-dependent 
coefficients a{t) and b{t) the original Emden equation can be reduced to a 
Lie system and then we can obtain a first integral which provides us with a 
i-dependent constant of the motion for the original system. 

In fact consider the system of first-order differential equations (fTOl) . This 
system describes all the systems of differential equations that can be ob- 
tained by means of the set of t-dependent transformations we got through 
the scheme S(WEmd,VEmd)- We recall that the t-dependent change of vari- 
able which we use to relate the Emden equation ([8]) with the latter system of 
differential equation is ([9]). As in previous papers about this topic we try to 
relate the initial system of differential equations to a Lie system determined 
by a t-dependent vector field of the form X'{t) = f{t)X and we suppose 
f{t) to be non- vanishing in the interval we study. So the system of differ- 
ential equations determining the integrals curves for this t-dependent vector 
field is a Lie system and we can use the theory of Lie systems to analyse its 
properties. 

As a first example we can consider that we just use the set of transfor- 
mations with 7(t) = 1 and a{t) =0. In this case system ^ is 



' dx' 
Id 




We fix j3{t) to be such that 



a{t) 



i.e., j3{t) is (proportional to) 
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Therefore we get 



^ = exp(^y a{t')dt'^ v', 

b{t) exp ^- J a{t')dt'^ x' 



In order to get the last system of differential equations to describe the integral 
curves for a f-dependent vector field, X'{t, x) — f{t)X{x), for a given function 
a{t) a necessary and sufficient condition is 



h{t) exp (^2 j a(t')dt'^ = K, 



with K being a real constant. Under this assumption the last system becomes 

dx' 



exp 



J\(t')dt'^ v', 



dt 

^ = exp( / a{t')dt')Kx"'. 



We introduce the t-reparametrisation 



it) = f^^p(^J a(t")dt"j dt' 



and the latter system becomes 




= V 



= Kx"", 



which admits a first integral 



1 



.X 



/n+l 



I = -v'^ - K- 

2 n + l 

In terms of the initial variables the corresponding t-dependent constant of 
the motion is 



/ = exp ( -2 / a(t')dt' ) ( -i' - b{t) 



n+l 
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which is similar to that found in 

Suppose that we restrict the transformations to the case a(t) = 0. In 
this case the system of first-order differential equations ffTOj) becomes 



dx' 
'dt 

d^ 



7(t) 7(t) 

m 



ait)- 



+ m ■ 



In order for this system of differential equations to determine the integral 
curves for a t-dependent vector field of the form X'{t) = f{t)X we need that 



lity 



C2i/(t) = ait) 



C22/(t) 



(33) 



From these relations, or more precisely from those of the first row, we get 
fit) as 

j^^^ ^ 1 lit) 1 m 



and therefore 



Cii 7(^) Ci2 7(^) 

Cl2 



We choose Cn = — 1 and Ci2 = 1 so that 

m = iit)- (34) 

In view of this and using the third and second relations from (!33|l we get 

ci2 7(t)" ^ ^ m 

and thus, as a consequence of the last differential equation becomes 



C21 7(^) ^ ^.^^^ _ 

ci2iit) at) 

and, as Ci2 = 1 and fixing C21 = 1, we obtain 



— log(77) = a(t). 
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which can be rewritten as 



Hence we have 



and in view of (134 



^^7'W=exp(|*a(t')dt'). 
7(t) = W2y" exp(^y a(t")dt"^ dt' 



f3{t) = ^ exp ( j a(t')dt'') . 

'2/*exp (^/*'a(t")dt") dt' ^ ^ 

So far we have used only three of the four relations we found. The fourth 
and second relations lead to the integrability condition: there exists a con- 
stant C22 = K such that 

m _ h{t)r{t) 
lit) m ■ 

Therefore, using the above expressions for 7(t) and [3{t), we get 

b{t)exp(^-2j a{t)dt'^ J exp^y a(t")dt"jj = K. (35) 

So under this assumption we have connected the initial Emden equation 
with the Lie system, 

d^ 

and then the method of characteristics shows that it admits the first integral 

2 n + 1 

In terms of the initial variables the corresponding constant of the motion is 



J = ( ia;2 _ i(^^n+i^ [\^t')dt') /*exp ( a(t")dt" | dt' 



^xxexp(^-j a(t')dt'^ (36) 
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and in this way we recover the result found in [5T]. If we now consider the 
particular case n = —3 we get that the integrability condition (1351) implies 
that there is a constant K such that 



b{t) 



exp (^-2 j a(t)dt'^ = K, 



and the corresponding t-dependent constant of the motion is then given by 
_^ !^iiZ.^-2 \ ( _2 / „,r^'^dt' W exn ( /* 



2 

which is equivalent to that found in [2T] . 



— -xiexp ( — / a(t')dt' ) , 



t 



8 A physical example. 

In order to motivate the physical applications of the theory of quasi-Lie 
schemes to Emden-Fowler equations, we here develop a detailed particular 
application of our methods to the Emden-Fowler equations of the form 

2 

x + -x + x" = 0, n 7^ 1. (37) 

t 

The case n = 1 is not dealt with in this work because quasi-Lie schemes are 
not needed to treat it. In this particular instance, the second-order differen- 
tial equation (I57|) is linear and the introduction of the new variable v = x 
transforms it into a nonautonomous linear system of first-order differential 
equations admitting a superposition rule. Thus fl37|) is a SODE Lie system 
[22] and the theory of quasi-Lie schemes reduces in this case to the usual 
theory of Lie systems. 

The equations (1571) initially appeared in the study of the stellar structure 
to analyse the configuration of spherical clouds of gas. These equations are 
also used in gas dynamics, fluid dynamics and more recently in relativistic 
mechanics, nuclear physics and chemical reacting systems [T7] . 

We must point out that the analysis of the solutions for equations (1371) is 
particularly important. Special interest have those non-oscillatory solutions 
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with positive zeros. These zeros are related to equihbrium states of a fluid in 
a spherical distribution of density and under mutual attraction of particles. 
Some results about the existence of these solutions can be found in |T71 [231 
[21]. Furthermore, those solutions x{t) with initial conditions a;(0) = 1 and 
x{0) = are specially important in astrophysics [25l [26] . In what follows, 
we analyse some transformation properties of the equations f[57l) and derive 
some new particular solutions. These solutions are interesting because they 
can provide information about stability of fluids in a spherical configuration 
or they can be used to check out methods to analyse particular solutions 
[231 [21]. 

It is known that fl37l) admits certain analytic solutions for the cases n = 
0, 1 and n = 5, e.g the case n = has a particular solution x(t) = — 
the case n = 1, the particular one x{t) = sint/t and the case n = 5 has 
the particular solution Xp{t) = (2t)~^/^, see |T7] and previous Sections. From 
the latter solution, we can recover another very well-known solution of the 
equation f[37|) . Suppose we want to get the particular solution x(t) for the 
equation (^7\} with x{0) = 1 and x(0) = 0. Thus, using this values at t = in 
the constant of the motion f[2Bl) obtained by means of the particular solution 
Xp{t) = (21)^^^"^, we get that such a constant is zero along the solution x{t) 
with the previous initial conditions and we have 

— — + At^x'^ + Ai^xx = 0, 

from which we obtain 

dx X I ^ / 4 , „\ 

— = — -1 + a/1 - -xH^ . 
dt 2ty y 3 J 

We make the change of variables x"^ = m and r = in order to simplify the 
above differential equation and to get 

dm m I r 4 ~\ 

— = — -1 + a/ 1 - -mV . 
dr y 3 J 

A new change of variables n = l/m transforms the above equation into 

dn 1 / n 

dr 2Ty V 3 ] 
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In order to remove the squared term in the previous differential equation, 
at least for a certain particular solution, we perform the change of variables 
n = u + r/3 and we obtain 



du I I ( T rZ ~r At\ 

— + - = — \u-\ \ + 2-u + — . 

rfr 3 2r \^ 3 V 3 3^ 3 ^ 

It is easy to check that the latter has a particular solution u = 1. Then 
inverting all previous changes of variables we get the well-known particular 
solution for (137|) with a;(0) = 1 and x{0) = 1 is given by 




We show that the knowledge of this particular solution allows us to obtain a 
way to get a family of exact solutions for the equation fl37j) with n = 5. In 
this way we illustrate how the knowledge of a particular solution can be used 
to obtain information about the general solution. From a physical point of 
view, this fact enables us to investigate by means of particular solutions the 
existence in the example of bounded oscillatory solutions the zeros of which 
determine equilibrium states of a fluid with spherical distribution of density 
and under mutual attraction of its particles [271 128] . 

Recall that the particular solution Xp{t) = (2^)"^/^ allowed us to trans- 
form the first-order system of differential equations related to (1371) to the one 
describing the integral curves for the t-dependent vector field 

= fit) {{x' + v')d,, - {v' + x")d,,) . 

Lie's theorem shows that any system determined by a t-dependent vector 
field on a manifold of the form X{t,p) = f{t)X{p), with p e A^, is a Lie 
system related to a onedimensional Lie algebra of vector fields. Moreover, 
the theory of Lie systems explains how to get a superposition rule for this Lie 
system [5] . Indeed, we have to consider m + 1 copies of the system. In this 
case, the solutions of the^ total system are integral curves for a t-dependent 
vector field of the form X{t) = f{t)X with 

m = fit) ^f2ix[^, + vi,)d,^^^ - + x[%)d..^^ . 
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The method developed in [5] explains that we have to consider m in such 
a way that vr^X be linearly independent with tt the canonical projection 
vr : (TR)"*^'""^-'^) — (TR)"™ onto the last nm variables. Then the first-integrals 
for X provide as the superposition rule. In our case, for m = 1 the vector 
field 7r*X is linearly independent and the first-integrals for X lead to the 
superposition rule. More explicitly, the vector field 

1 

a=0 

admits locally three first-integrals allowing us to obtain the superposition 
rule. The method of characteristics provides us with two of such first-integrals 



1 'g 1 '2 / / 



0,1. 



In order to obtain an analytic expression for a third (partial) first-integral, 
we restrict ourselves to the submanifold 

C = {{x',v') e {TRf\Io{x',v') = h{x',v') = 0}, 

with {x',v') = (x|qj, f ^Q-), f ^-1^^). As is shown below, this case is enough 
general to obtain new exact solutions for (P7|) . Under the assumed restric- 
tions, a third first-integral K on C, that is XK\c = 0, is 



K 




(38) 



From the latter expression we can work out the value of x'^q^ in terms of x'^^-^ 
and the constant K 



X 



(0) 



\ 



QKx'l^^ 



„/4 

3 



1 



llii 

3 



12K^ + (1 - K2) V(^)^ 



and taking into account that x'^-s^ = y'^tXi, with i = 0, 1, we get what we call 
a partial t-dependent superposition rule for the equations ( |371) : 



Xq 



\ 



6Kxf 1 - V 1 



3 



1 



3 



12^2 + (1 - 7^2)24^23,4 
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It can be seen that this t-dependent partial superposition rule does not gen- 
erate all solutions for the Emden- Fowler equation (1371) . The reason is that in 
order to get an algebraic expression for such a partial superposition we had to 
restrict ourselves to some kind of solutions. The name of partial is related to 
the partial superposition rule concept [3]. These partial superposition rules 
do not enable us to obtain the general solution, but a family of solutions. 

If we take into account that Xi{t) = (1 + t^/3)~^/^ is a particular solution 
for the equation (1371) with n = 5, we can use such a solution in our t- 
dependent partial superposition rule to get the K parametrised family of 
solutions 



Xo{t) 




3 + t2| + K2(3 + t2 + I - 3 + t2|)) 



(3A^2 + ^2)(3 + ^2^2) 



In particular, for = 1 we recover Xoit) = (1 + and for K = we 
get Xo(t) = 0. We can see that for any K we get that Xq{0) = 0. 

This t-dependent partial superposition rule can be used to analyse the 
solutions of the equation fl37|) . In particular, it can be seen that in the limit 
t — > cx) all solutions xo(t) tend to zero, i.e. lim^^oo a;o(t) = 0. 

Finally, we proved in Section [6] that the system ( l30i) related to an Emden- 
Fowler equation of the form ( l29l) can be transformed by means of a t- 
dependent change of variables Q) with a(t) = j(t) and P and 7 two t- 
dependent functions holding the differential equations 

7 = -q{th - Pm and - p{t) - ^ " §^ = 0, 

lit) Pit) 

into the system 

dx' _(3{t) , 

d& ^ r{t)r{t) 
dt Pit) ■ 

For the case of Emden- Fowler equations fl37|) we have that p{t) = 2/t, q{t) = 
and r(t) = 1. Thus, we get that 

7=A and -i-m_m=o 
^ t lit) Pit) 
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Hence, jit) = P{t) = t ^ and we obtain 




v' 



Therefore, we can see that for any solution {x' (t) , v' (t)) of the above system, 
the function x'{t) satisfies the second-order differential equation 



a2 / 



and the change of variable x = t~^x' transforms the equation (I57|) into the 
above differential equation. Even whether this result can be found in [T7], 
we here provide a new approach to it. 



9 Conclusions and Outlook. 

We have applied the theory of quasi-Lie schemes to treat Emden equa- 
tions. We have shown that the knowledge of a particular solution satisfying 
a certain condition enables us to transform a given Emden equation into a 
Lie system. As a consequence we get a family of Emden equations which are 
quasi-Lie systems with respect to the quasi-Lie scheme we have used. This 
fact provides a method to obtain t-dependent constants of the motion. We 
have offered a method to construct Emden equations with a previously fixed 
particular solution related to a t-dependent constant of the motion. We have 
also explained the Kummer-Liouville transformation from the point of view 
of the theory of Lie schemes. Next, we have applied a method used in previ- 
ous papers about quasi-Lie schemes to obtain t-dependent constants of the 
motion for certain families of Emden equations. Finally, in order to motivate 
the physical applications of our methods, we have analysed certain cases of 
Emden-Fowler equations. For such cases, some properties of transformation 
have been derived. Additionally, for an interesting particular Emden equa- 
tion, we have recovered a particular solution and through such a solution, 
and by means of a t-dependent partial superposition rule for this Emden 
equation, a family of, as far as we know, new solutions has been derived. 

As it was stated in |6J, the field of apphcations of quasi-Lie schemes is 
very broad. Here we have presented a new application of this concept and we 
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hope that more results can be recovered or obtained by means of our methods, 
e.g. new t-dependent partial superposition rules to obtain particular exact 
solutions from particular ones. 
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